■  ESavntHS8«£6  _  j_:% 

P  ESP  RECORD  COPY 

W-'iy  ’  '  J  : 

■  SCIENTIFIC  ft  TECHNICAL  ISFORMATIOH  0IV8WH 


\ 


(ESTI),  BUILDING  1211 


ESD  ACCESSION  USE 

ESTI  Can  No.  H  60861  ’ 

Copy  No,  _L^L*_ZT£ 


* 


Technical  Note 


1967-16 

Rev.  1 


J.  U.  Beusch 


Stable  Equilibrium  Orientation 
of  a  Spacecraft  Which  Contains 
a  Momentum  Wheel 


10  May  1968 


Prepared  under  Electronic  Systems  Division  Contract  AF  19(628K5167  by 

Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexington,  Massachusetts 


The  work  reported  in  this  document  was  performed  at  Lincoln  Laboratory, 
a  center  for  research  operated  by  Massachusetts  Institute  of  Technology, 
with  the  support  of  the  U.S.  Air  Force  under  Contract  AF  19(628)-5167. 

This  report  . may  be  reproduced  to  satisfy  needs  of  U.S.  Government  agencies. 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 


LINCOLN  LABORATORY 

STABLE  EQUILIBRIUM  ORIENTATION  OF  A  SPACECRAFT 
WHICH  CONTAINS  A  MOMENTUM  WHEEL 

J.  U.  BEUSCH 
Group  67 


TECHNICAL  NOTE  1967-16,  REV.  1 

10  MAY  1968 


LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


The  orientation  of  a  spacecraft  when  it  is  in  a  stable  equilibrium 
state  is  studied.  The  spacecraft  contains  a  single  momentum  wheel  which 
stores  angular  momentum.  The  axis  of  the  wheel  is  not  necessarily  parallel 
to  a  principal  axis  of  the  spacecraft.  For  an  arbitrary  but  fixed  speed  of 
the  wheel  relative  to  the  spacecraft  body,  it  is  shown  that  there  may  be  one, 
two,  or  three  stable  equilibrium  points.  At  each  of  these  points,  the  body 
may  spin  about  an  axis  which  is  fixed  in  both  spacecraft  and  inertial  coor¬ 
dinates.  The  orientation  of  this  axis  in  spacecraft  coordinates  can  be 
determined  from  expressions  in  this  report.  A  threshold  can  be  deter¬ 
mined  such  that,  if  the  speed  of  the  wheel  relative  to  the  spacecraft  is 
larger  than  this  threshold,  there  is  only  one  stable  equilibrium  point.  This 
information  can  be  used  to  determine  the  eventual  orientation  of  the  space¬ 
craft  provided  it  is  designed  such  that,  if  it  becomes  seriously  misoriented, 
the  momentum  wheel  drive  motor  automatically  holds  the  wheel  speed  to  a 
predetermined  value. 
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STABLE  EQUILIBRIUM  ORIENTATION  OF  A  SPACECRAFT 


WHICH  CONTAINS  A  MOMENTUM  WHEEL 


1.  Introduction 

The  orientation  of  a  spacecraft  when  it  is  in  a  stable  equilibrium  state 
is  studied  in  this  report.  The  spacecraft  is  assumed  to  contain  a  single 
momentum  wheel  with  its  spin  axis  in  an  arbitrary  direction  and  a  damper  for 
dissipation  of  energy.  Two  situations  are  already  well  understood. 

(a)  The  wheel  is  spinning  at  such  a  speed  that  the  magnitude  of 
its  angular  momentum  equals  that  of  the  entire  spacecraft.  In  its  unique  stable 
equilibrium  state,  the  spacecraft  will  orient  itself  such  that  the  momentum 
vector  of  the  wheel  coincides  with  the  momentum  vector  of  the  entire  space¬ 
craft  which  is  fixed  in  inertial  coordinates. 

(b)  The  wheel  has  stopped  spinning  relative  to  the  spacecraft  body. 
There  are  two  stable  equilibrium  states.  In  each,  the  axis  of  maximum  moment 
of  inertia  of  the  spacecraft  is  oriented  in  one  of  two  possible  directions  along 
the  momentum  vector  and  the  spacecraft  is  spinning  about  this  axis. 

In  each  of  these  two  situations,  when  the  spacecraft  is  not  in  a  stable 
equilibrium  state,  the  damper  dissipates  energy  and  drives  the  spacecraft 
toward  a  stable  equilibrium  state.  *  In  this  stable  equilibrium  state,  the  damper 
dissipates  no  energy  since  the  spacecraft  is  either  stationary  or  spinning  about 
a  fixed  axis. 

Intermediate  situations  between  these  two  extreme  situations  are 
analyzed  in  this  report.  The  momentum  wheel  speed  relative  to  the  body  is 
assumed  to  be  held  constant  by  a  drive  motor  and  the  orientation  of  the  space¬ 
craft  when  it  is  in  a  stable  equilibrium  state  is  investigated.  It  will  be  shown 
that  for  an  arbitrary  but  fixed  wheel  speed  in  any  stable  equilibrium  state  the 

*  In  theory,  no  dissipation  takes  place  when  the  spacecraft  is  in  an  unstable 
equilibrium  state,  but  this  case  is  ignored  in  this  report  because  in  practice 
a  disturbance  torque  will  perturb  the  spacecraft  from  an  unstable  equilibrium 
state. 
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body  will  rotate  about  some  axis  which  orients  itself  along  the  spacecraft 
momentum  vector.  The  damper  'dissipates  no  energy  while  the  spacecraft 
spins  about  this  axis.  This  axis  is  not  necessarily  the  wheel  axis  or  a  princi¬ 
pal  axis  of  the  spacecraft.  If  the  wheel  speed  is  larger  than  some  threshold, 
it  will  be  shown  that  there  is  only  one  stable  equilibrium  state  and  the  orienta¬ 
tion  of  the  spacecraft  when  in  this  state  can  be  determined.  If  the  spacecraft 
is  designed  so  that  when  it  becomes  misoriented  the  momentum  wheel  is  auto¬ 
matically  driven  at  a  predetermined  speed,  the  methods  of  this  report  can  be 
used  to  determine  the  equilibrium  orientation  of  the  spacecraft. 

The  precise  evaluation  of  the  expressions  of  this  report  requires  a 
knowledge  of  h,  the  magnitude  of  the  spacecraft  angular  momentum.  For 
example,  the  threshold  varies  linearly  with  h.  However,  in  many  cases,  the 
direction  of  the  axis  in  spacecraft  coordinates  about  which  the  body  spins  is 
insensitive  to  fairly  large  variations  in  h  as  long  as  the  wheel  speed  remains 
above  the  threshold.  Therefore,  the  results  of  this  report  are  useful  when  h 
is  known  only  approximately. 

When  two  of  the  principal  moments  of  inertia  are  equal,  the  results  of 
this  report  require  some  special  interpretation  since  some  of  the  expressions 
become  infinite.  The  best  way  to  handle  this  situation  is  to  let  the  moments  of 
inertia  differ  by  a  small  quantity  and  observe  what  happens  when  this  quantity 
approaches  zero.  In  practice,  no  two  moments  of  inertia  are  exactly  equal, 
so  the  results  for  equal  moments  of  inertia  are  not  emphasized  in  this  report. 

For  the  special  case  where  the  spacecraft  body  is  symmetric  about  the 
wheel  axis,  the  threshold  value  of  the  wheel  speed  is  determined  in  Ref.  2 
using  a  somewhat  different  approach. 

2.  Formulation  of  the  Problem 

Consider  a  spacecraft  which  is  rigid  except  that  it  contains  a  momentum 
wheel  that  is  free  to  rotate  about  a  shaft  parallel  to  an  axis  Z  in  spacecraft 
coordinates  which  passes  through  the  spacecraft  center  of  mass.  Define  the 
direction  of  positive  Z  such  that  a  vector  representing  angular  velocity  of  the 
wheel  with  respect  to  the  rest  of  the  spacecraft  is  in  the  direction  of  positive  Z. 
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The  momentum  wheel  is  assumed  to  have  symmetry  about  its  spin  axis.  Two 

additional  coordinates  X  and  Y  are  chosen  so  that  XYZ  forms  a  right-hand 

orthogonal  coordinate  system  centered  at  the  spacecraft  center  of  mass  and 

fixed  to  the  rigid  part  of  the  spacecraft  which  does  not  include  the  wheel.  The 

wheel  can  be  considered  to  be  composed  of  a  first  part  with  zero  inertia  about 

axes  X  and  Y  and  inertia  C  about  Z  and  a  second  part  with  zero  inertia  about 

axis  Z  and  non-zero  inertia  about  axes  X  and  Y.  The  first  part  is  free  to 

rotate  about  Z  relative  to  the  rest  of  the  spacecraft.  The  second  part  also 

rotates  with  the  first  but  for  purposes  of  obtaining  an  expression  for  the 

kinetic  energy  it  is  convenient  to  assume  that  the  second  part  is  rigidly  attached 

to  the  spacecraft.  This  represents  no  restriction  since  the  kinetic  energy  is 

the  same  regardless  of  the  rotational  speed  about  Z  of  the  second  part.  All 

parts  of  the  spacecraft  other  than  the  first  part  of  the  wheel  will  be  referred 

to  as  the  spacecraft  body  and  assumed  to  be  rigid. 

The  spacecraft  body  has  three  principal  axes.  Denote  these  axes  x, 

y,  and  z,  which  form  a  right-hand  orthogonal  coordinate  system  such  that  the 

positive  Z  axis  is  in  the  octant  x  ^  0,  y  £  0,  and  z  £  0,  and  such  that  I  ^  I 

1  x  y 

and  I  ^  I  where  1,1,  and  I  are  the  principal  moments  of  inertia  of  the 
x  z  x  y  z 

spacecraft  body.  See  Fig.  1. 

Since  the  wheel  has  zero  inertia  about  X  and  Y  its  angular  momentum 
is  a  vector  go  C  along  Z  where  go  is  the  component  of  angular  velocity  of  the 
wheel  along  Z.  The  spacecraft  has  a  total  angular  momentum  h  in  some 
arbitrary  direction.  The  vector  h  can  be  broken  into  two  components,  a 
vector  Cgo  along  Z  and  a  vector  h  —  Cgo.  Let  the  components  of  Cgo  along  x, 
y,  and  z  have  magnitude  a^  Cgo,  a^  Cgo,  and  a^  Cgo  where  go  =  |go|  and  let  the 
components  of  h  have  magnitude  b  h,  b  h,  and  b  h  where  h  =  jh  j  .  Then 

zzz  z  z  z  x  y  z 

a  ^  +  a  +a  =1,  b  +  b  +  b  »  1,  and,  because  of  the  way  x,  y,  and  z 
x  y  z  x  y  z  7  7 

were  defined,  a  £  0,  a  £  0,  and  a  ^  0. 

x  y  z 

3.  Points  of  Minimum  Kinetic  Energy 

In  this  section,  we  consider  a  problem  where  a,  a,  a,  I,  I,  I, 

x  y  z  x  y  z 

C,  h,  and  go  are  fixed  and  examine  the  minima  of  the  function 


ZT(b 


x 


bz)A2  =  (bx-pax)2/lx  +  (b-cay)'i/ly  +  (bz-Oaz)2/lz+p2/C 


(1) 
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CO  Q 

where  T  is  the  kinetic  energy  of  the  spacecraft  and  p=  (— r—  )•  Equation  1  follows 

2  “ 

from  the  fact  that  the  energy  of  the  wheel  is  (go  C)  /ZC  and  the  energy  of  the  body 
is 


(hb  -a  wC)Z  (hb  -a  u  C)Z  (hb  -a  wC)Z 

x  x  +  y  y  .  z  z 

21  21  21 

x  y  z 

2  2  2 

For  each  point  on  the  unit  sphere,  b  *fbr  lb  =1,  the  value  of  T  can 

x  y  z 

be  computed  using  Eq.  1.  The  assumption  that  go  is  fixed  implies  that  a  motor 
drives  the  wheel  such  that  the  component  of  its  angular  velocity  along  Z  is  con¬ 
stant.  Assuming  that  energy  is  dissipated  in  the  body  by  the  damper,  regardless 

of  the  initial  values  of  b  ,  b  ,  and  b  (initial  orientation  of  h  in  the  spacecraft 

x  y  z  — 

coordinate  system),  the  spacecraft  will  eventually  orient  itself  so  that  b^,  b  , 

and  b  take  on  the  values  of  the  coordinates  of  a  minimum  of  T.  In  this  stable 
z 

equilibrium  state,  the  direction  of  h  is  fixed  in  both  inertial  coordinates  and  in 
the  spacecraft  coordinates.  Therefore,  the  spacecraft  body  is  either  stationary 
with  respect  to  an  inertial  frame  of  reference  or  rotating  about  h.  In  particular, 
its  component  of  angular  velocity  along  Z,  Q^9  is  constant.  In  this  report  the 
terms  minimum  and  stable  equilibrium  point  are  used  interchangeably.  At  a 
minimum,  the  damper  dissipates  no  energy  since  the  body  is  either  stationary 
or  spinning  about  a  fixed  axis. 

In  the  above  formulation  it  is  assumed  that  go,  the  component  along  Z  of 

the  angular  velocity  of  the  wheel  relative  to  inertial  coordinates,  is  constant. 

In  practice  the  determination  of  stable  equilibrium  points  when  the  speed  of  the 

wheel  relative  to  the  body,  go  =  go  —  0  ,  is  constant  is  of  most  interest.  However, 

the  solutions  to  the  two  problems  are  the  same  as  long  as  the  variables  are 

interpreted  properly.  This  is  true  because  at  equilibrium  is  constant  as 

stated  in  the  preceding  paragraph.  In  many  practical  cases  Q  is  small  compared 

to  go  so  that  go  «  go  . 

s 

There  are  several  situations  in  which  knowledge  of  the  components  of  h, 
or  equivalently  the  direction  of  the  axis  about  which  the  body  spins,  in  spacecraft 
coordinates  is  valuable.  Two  situations  are: 

(a)  The  spacecraft  loses  power  for  a  while  (e.  g.  ,  while  in  the 
earth’s  shadow)  so  that  the  momentum  wheel  stops  spinning  relative  to  the  body. 
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Power  is  regained  and  the  momentum  wheel  is  automatically  accelerated  to  a 
predetermined  speed  o>s  relative  to  the  body.  In  many  practical  cases,  the 
magnitude  of  the  angular  momentum,  h,  will  be  much  larger  than  the  magnitude 
of  angular  momentum  changes  caused  by  disturbance  torques  during  the  power 
outage  so  that  the  direction  and  magnitude  of  h  remain  approximately  fixed  in 
inertial  space.  Knowledge  of  the  components  of  h  in  the  spacecraft  coordinate 
system  would  be  extremely  helpful  in  bringing  about  complete  recovery  of  the 
spacecraft  from  misorientation  suffered  while  the  momentum  wheel  was  without 
power. 

(b)  Due  to  faulty  ejection  of  the  spacecraft  from  the  launch  vehicle 
or  large  disturbance  torques  caused  by  failure  of  gas  jets  in  the  open  position, 
the  magnitude  and  orientation  of  h  in  inertial  coordinates  is  unknown  but  the 
wheel  is  automatically  held  by  the  motor  to  a  predetermined  speed  oo  relative  to 
the  body.  Knowledge  of  the  orientation  of  h  in  spacecraft  coordinates  would  be 
helpful  as  a  first  step  in  determining  the  orientation  of  h  in  inertial  space  and 
eventual  spacecraft  recovery. 


It  is  shown  in  the  appendix  that  the  function  T  (b^,  b^,  b^)  always  has 

one  and  only  one  minimum  in  the  region  b  ^  0,  b  ^  0,  b  ^  0.  T  may  have 

x  y  z 

zero,  one,  or  two  minima  in  the  region  b  <0,  b  ^  0,  b  ^  0  and  there  are 

x  y  z 

no  minima  in  any  other  regions.  Two  sufficient  conditions  that  T  have  a 
unique  minimum  are 


<i£>2  > 


(I  -I  ) 

x  y 


(a  I  )2/3  +  (a  I  )2/3  3 
x  y  y  x  j 


and 


(^)2  > 


(I  -I  ) 

x  z 


L( 


a  I  ) 

X  z 


+  (a  I  )2/3 

Z  X 


U) 


(3) 


If  either  Eq.  2  or  3  is  satisfied,  T  will  have  one  and  only  one  minimum.  Let 

t  =  C/h  times  the  positive  square  root  of  the  smaller  of  the  two  expressions  on 

the  right-hand  side  of  Eqs.  2  and  3.  The  coordinates  of  this  minimum  satisfy  the 

relation  given  in  Table  1  for  various  ranges  of  the  parameters  a  ,a  ,  and  a  . 

x  y  z 
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Table  1 


Conditions 

on  a  ,  a  ,  a 
x  y  z 

Coordinates  of  Unique  Minimum  Satisfy 

a  =1, 

a 

_ 

a  =  0* 

b 

1, 

b  =  b  =0 

X 

y 

z 

X 

y  z 

a  =1, 

a 

— 

a  =  0* 

b 

— 

1, 

b  =  b  =  0 

y 

X 

z 

y 

X  z 

a  =1, 

a 

= 

a  =  0* 

b 

— 

1, 

b  =  b  =  0 

z 

X 

y 

z 

x  y 

a  =0, 

a 

> 

0,  a  >  0* 

b 

- 

0, 

b  =  pa  I  b  /[  pa  I  +  (I  -I  )b  ]  , 

X 

y 

z 

X 

y  yzz  zy  zyz 

b 

= 

+ 

1  — b  2 

y 

<v 

z 

a  =0, 

a 

> 

0,  a  >0 

b 

- 

0, 

b  =  oa  I  b  /[pa  I  +  (I  -I  )b  ]  , 

y 

X 

z 

y 

z  zxx  xz  xzx 

b 

- 

+  / 

1-b  2 

X 

-v 

z 

a  =0, 

a 

> 

0,  a  >0 

b 

- 

0, 

b  =  Pa  I  b  /[pa  I  +  (I  -I  )b  ]  , 

z 

X 

y 

z 

y  y  x  x'  x  y  x  y  xJ 

b 

= 

+  7 

1  — b  2 

X 

V 

y 

a  >0, 

a 

> 

0,  a  >0 

b 

= 

pa 

I  b  /[pa  I  +  (I  -I  )b  ] 

X 

y 

z 

z 

z 

XX  X  Z  X  Z  X 

b 

= 

pa 

I  b  /[pa  I  +(I  -I  )b  1 

y 

y 

xx  x  y  x  y  x 

b 

_ 

+ ! 

l_b  2  — b  2 

X 

y  z 

'  ^  these  four  cases,  a  necessary  condition  for  T  to  have  a  unique  minimum 
is  that  either  Eq.  2  or  3  with  >  replaced  by  2:  be  satisfied. 
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Given  the  values  of  the  spacecraft  parameters,  a  threshold  t  of  wheel 
speed  can  be  determined  from  Eqs.  2  and  3  such  that  if  cj  >  t  the  vector  h 
has  only  one  stable  equilibrium  point  in  the  spacecraft  coordinate  system.  In 
this  equilibrium  position,  the  components  of  h  can  be  determined  from  the 
relations  of  Table  1.  If  cj  <  t,  there  may  be  as  many  as  three  equilibrium 
positions.  In  each  equilibrium  position,  the  spacecraft  body  may  rotate  about 
h.  These  positions  may  be  studied  using  the  expressions  of  the  appendix.  It 
is  desirable  that  the  spacecraft  be  designed  to  keep  oa  >  t  so  that  the  additional 
equilibrium  positions  will  not  exist. 

4.  Example:  Wheel  Axis  Along  a  Principal  Axis 

Consider  a  spacecraft  composed  of  two  rigid  bodies.  See  Fig.  2.  The 
first  body  contains  a  momentum  wheel  with  its  spin  axis  parallel  to  Z.  The 
second  body  rotates  very  slowly  (e.  g.  ,  once  per  day  for  spacecraft  in  synchro¬ 
nous  equatorial  orbit)  about  Z  with  respect  to  the  first  and  Z  is  a  principal 
axis  of  both  bodies.  This  slow  rotation  is  caused  by  a  drive  which  for  long 
periods  of  time  (e.  g.  ,  twenty-four  hours)  does  not  require  any  sensor  input 
from  outside  the  spacecraft  so  that  during  any  particular  short  time  interval 
of  interest  the  spacecraft  can  be  considered  to  be  a  rigid  body  with  one  prin¬ 
cipal  axis  along  axis  Z.  Assume  that  the  moments  of  inertia  normalized  by 
dividing  by  the  constant  1^  are 

I  /lz  =  5.  330  +  0.  038  cos  20, 

I  /i  =  5.  160  -  0.  038  cos  28, 

1  z 

and 

I^y/^Z  =  °‘  38  SiT1  26 

where  0  is  the  angle  of  body  1  relative  to  body  2. 
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Consider  the  satellite  position  9  =  0.  In  this  case,  the  X,  Y,  and  Z 
axes  are  the  principal  axes, 


and 


I  /i  =  5.  368, 

X  z 


I  /i  =5.  122. 

y  z 


The  case  described  here  is  =  1,  ax  =  ay  =  Equations  Z  and  3  imply  that, 

if  0  =  C^-)  >  1— (I  /i  )  =  0.  814,  a  unique  minimum  occurs  at  b  =  1,  b  = 
n  z  x  z  y 

b^  =  0.  Therefore,  if  the  speed  of  the  wheel  is  always  kept  high  enough  so 
that  81.  4$  of  the  angular  momentum  of  the  spacecraft  is  due  to  the  wheel 
spinning,  only  one  stable  equilibrium  position  exists.  In  this  position,  the 
body  may  be  spinning  around  the  Z  axis. 

If  o  <  0«  814,  the  results  of  the  appendix  show  that  there  may  be  two 
equilibrium  positions  at 

b  =  0, 

y 


bz  =  DIx/(Ix-Iz)  =  p/0.  814, 


and 


bx  =  ±.p-(p/0.  814)2. 

Consider  the  satellite  position  0  =  tt/Z.  Equations  Z  and  3  imply  that, 

Q 

if  o  =  (-1 — )  >  0.  811,  a  unique  minimum  occurs  at  b  =  1,  b  =  b  =0.  At  this 
M  h  ^  z  y  x 

stable  equilibrium  position  the  body  may  spin  about  the  Z  axis. 

5.  Example:  Wheel  Axis  Which  may  not  be  Along  a  Principal  Axis 


Consider  a  spacecraft  composed  of  two  rigid  bodies.  See  Fig.  3.  The 
first  body  contains  a  momentum  wheel  with  its  spin  axis  parallel  to  Z.  The 
second  body  rotates  with  respect  to  the  first  very  slowly  about  Y  and  Y  is  a 
principal  axis  of  both  bodies.  With  respect  to  this  rotation,  the  situation  is 
the  same  as  that  of  the  previous  example  so  that  for  stability  analysis  the 
spacecraft  can  be  considered  to  be  a  single  rigid  body  which  contains  a 
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momentum  wheel.  Assume  that  the  moments  of  inertia,  normalized  by 
dividing  by  the  constant  I are 

Ix/ly  =  7.  4  +  0.  089  cos  26, 

Iz/ly  =  6.  89  -  0.  089  cos  20, 

and 

lXZ^lY  =  °'  089  Sin  20 


where  0  is  the  angle  between  the  two  bodies  which  is  considered  to  be  a  con¬ 
stant  at  any  particular  time.  The  coordinates  x,  y,  and  z  are  chosen  so  that 
z  is  along  Y  and  x  is  in  the  XZ  plane  displaced  an  angle  cp  from  X.  The 

condition  I  =0  determines  cp.  By  substituting  coordinate  transformations 
xy  7  & 

into  the  definition  for  the  moments  of  inertia,  we  obtain  the  relations 


and 


ly  Xy(IZ  +IX)  +  (IX  lZ  ZXZ  ^ 


Xx  (IX!Z  ZXZ 


sm  cp  = 


ly-1 

_ y_ 

I  -  I 

X  y 


0, 


which  can  be  solved  for  I  /i  ,  I  /i  ,  and  cp. 

y  z  x  z 

Consider  the  case  0  =  rr/4.  Then 


I  /i  =6.  875 
Y  z 

I  /i  =  7.  415 
X  z 

and 

cp  =  9.  6°. 


a  =  0.  167,  a  =  0.  985. 
x  y 


Equations  Z  and 


Therefore,  a  =0, 

C  ^ 

p  =  (-^-)  >  0.  05,  there  is  only  one  point  of  minimum  energy. 


VTP 

b  =0  and 
x 


b  and  b  satisfy 
y  x  J 

b  =  6.  35  ob  /(p  +  0.  47b  ) 
y  M  x  x 


and 


b 

x 


3  imply  that,  if 
At  this  point, 


(4) 

(5) 
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This  follows  from  Table  1.  Therefore,  if  the  wheel  speed  is  kept  high  enough 
so  that  more  than  5#  of  the  angular  momentum  of  the  spacecraft  is  due  to  the 
wheel  spinning,  the  momentum  vector  has  only  one  stable  equilibrium  point  in 
spacecraft  coordinates.  In  this  equilibrium  position,  the  vector  h  is  not  parallel 
to  the  wheel  axis.  The  body  may  spin  about  an  axis  along  h. 

If  p  <  0.  051,  there  may  be  one,  two,  or  three  equilibrium  points 

depending  upon  whether  Eqs.  4  and  5  have  one,  two,  or  three  solutions  in 

the  region  b  =0,  1  ^  b  ^  0,  l^b  ^  -1.  This  follows  from  the  appendix. 

»  z  y  x 

Consider  the  case  0  =  0.  Then 
I  /i  =  7.  489, 

X  z 

I  /i  =6.  801, 
y  z 

and  cp  =  0°  so  that  a  =1  and  a^  =  a^  =  0.  If  p  >  0.  092,  there  is  only  one 
point  of  minimum  energy  located  at  b  =  1,  b^  =  =  0*  In  this  stable  equi¬ 

librium  state,  the  body  may  spin  about  the  Z  axis. 

If  p  <  0.  092,  there  are  two  stable  equilibrium  points  at  b^  =  0,  b  = 
p/o.  092,  and  =  ±  ^  1  —  (D /0.  092)^. 

Consider  the  case  0  =  tt/2.  Then 

I  /i  =  7.  31 1, 

X  z 


and  cd  =  0°  so  that  a  =1  and  a  =  a  =0.  If  n  >  0.  0454,  there  is  onlv  one 

y  x  z  7 

point  of  minimum  energy  located  at  b  =  1,  b  =  b  =0.  In  this  stable  equi- 

07  y  x  z  ^ 

librium  state,  the  body  may  spin  about  the  Z  axis.  If  p  <  0.  0454,  there  are 
two  stable  equilibrium  points  at  b^  =0,  b^  =  p/0.  0454  and  -b^  =  ± ^1  —  (p/0. 0454)^. 

6.  Conclusions 

For  t  <  co  <  h/C  where  t  is  some  threshold,  u>  is  the  momentum  wheel 
speed,  C  is  the  wheel  inertia,  and  h  is  the  magnitude  of  the  spacecraft  angu¬ 
lar  momentum,  there  is  only  one  equilibrium  state  of  the  spacecraft.  If  the 
momentum  wheel  axis  is  not  parallel  to  a  principal  axis  of  the  spacecraft,  in 
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the  equilibrium  state  the  spacecraft  body  spins  about  an  axis  which  is  not 
parallel  to  either  the  wheel  axis  or  a  principal  axis  when  t  <  co  <  h/C.  If  the 
wheel  axis  is  parallel  to  a  principal  axis,  in  the  equilibrium  state  the  body 
spins  about  an  axis  which  is  parallel  to  the  wheel  axis  when  t  <  go  <  h/C.  The 
value  of  t  and  the  orientation  of  the  spin  axis  of  the  spacecraft  body  in  the 
equilibrium  state  can  be  determined  using  the  expressions  derived  in  this 
report. 
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APPENDIX 


1.  All  minima  of  T(b^,  b^,  b^)  satisfy  the  following  expressions.  If 


pa  =0, 
x 


b  =  0 
x 


pa  I 

o  <  b  £  y  x 


y  I  —I 
’  x  y 


(A-  1 ) 


pa  I 

u  _  y  x 

y  "  (I  -I  ) 

1  x  y 


-  l  «;  b  <l. 
x 


(A-2) 


If  pa  =  0, 

y 


b  =  0 

y 


b  =  0 

y 


0  <  b  <1 
x 


-pa  I 

.  ,  x  y 

-  1  b  ^  -j - D- 

x  I  —  I 
x  y 


(A-  3) 

(A  -4) 


If  pa  >0  and  pa  >  0, 
x  y 


pa  I  b 
y  x  x 

V  =  pa  I  +  (I  -I  )b 
’  x  y  x  y  x 


0  ^  b  <,  \ 
x 


(A-  5) 


pa  I  b  /a  I  \l/3 

b  =  -  y  V  T  ,u  _l^b^-b(-V 

y  pa I  +  (I  —  I  )b  x  y\  a  ^  1^ 


x  y  x  y  x 


y 


-pa  I 

and  b  ^  -j — X  —  . 
x  I  —  I 
x  y 


(A-6) 


These  expressions  denote  surfaces  in  b  ,  b  ,  b  space.  On  the  b  ,b  plane 
r  x  y  z  ^  xy 

they  denote  curves  which  are  shown  in  Figs.  A-l,  A-2,  A-3,  and  A-4.  To 
derive  these  expressions,  let 


where 


T(b  ,  b  )  =  T./h2  +  T,/h2 
x  y  1  2 


T.  =  (b  -pa  )2/l  +  (b  -Da  )2/l 

l  x  w  x  '  x  y  y  y 
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and 


=  W1”11*  -  V  “  pa^  /xz  +p  /c* 


Define  r  and  0  such  that  b  =  r  cos  0  and  b  =  r  sin  0,  and  consider 

x  y 

circles  of  constant  r  on  the  b^,  b^  plane.  A  necessary  condition  for  a  point 

at  radius  r  to  be  a  minimum  of  T(b  ,  b  )  is  that  on  a  circle  of  constant  r 

x  y 

this  point  be  a  minimum  of  T(0),  a  function  of  a  single  variable  0.  Since  T  } 
is  independent  of  0,  a  point  is  a  minimum  of  T(9)  if  and  only  if  it  is  a  mini¬ 
mum  of  T^(0).  The  conditions 


i  dTi 


Z  d0 


b  (b  —  pa  ) 

x  y  y 


b  (b  —pa  ) 
y  x  x 


=  0 


and 


i  d'xi 


2  flC2 


-b  (b  —pa  )  +  (b  ) 

y  y  y  x 

i 


b  (b  -pa  )  —  (b  ) 
xxx  y 

I 


;>  0 


which  must  be  satisfied  by  all  minima  for  r  >  0  yield  the  desired  expressions. 
The  point  r  =  0  must  be  checked  separately. 


Z.  All  minima  of  T(b  ,  b  ,b  )  satisfy  the  following  expressions.  If 

x  y  z 


pa  =  0, 
x 


b  =  0 
x 


oa  I 

o  <  b  <;  - — 

z  I  —  I 
X  z 


(A-7) 


pa  I 
z  x 

(I  -I  ) 

X  z 


-  1  <  b  <1 
x 


(A-  8) 


If  pa  =0, 
z 


b  =  0 
z 


0  <  b  <1 
x 


(A-9) 


b  =  0 
z 


-oa  I 
x  z 

-  1  <  b  <  -J - y - 

x  1  —  I 
X  z 


(A-  1  0) 


If  pa  >0  and  oa  >0, 
x  z 
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(A  -  1  1 ) 


b 

z 


pa  I  b 

Z  X  X 

pa  I  +  (I  -I  )b 

X  Z  X  Z  X 


0  <  b  <  1 

X 


b 

z 


pa  I  b 

z  x  x 

pa  I  +  (I  -I  )b 

X  Z  X  Z  X 


-  1  <;  b  £  -b 

X  z 


1/3 


-pa  I 

,  .  x  z 

and  b  <  -= - = — 

x  I  —  I 
x  z 


(A- 1 2) 


These  expressions  are  identical  to  those  of  Section  1  except  that  y  is  replaced 
by  z.  They  are  derived  in  the  same  way  as  the  expression  of  Section  1. 


3.  Combining  the  results  of  Sections  1  and  Z  yields  curves  in  three- 

dimensional  space  on  which  all  minima  of  T  must  lie.  These  curves  are 

intersections  of  the  surfaces  of  Sections  1  and  Z  and  are  shown  in  Figs.  A-5 

through  A- 10  for  six  cases  that  occur.  A  seventh  case,  a  =  0,  a  >0,  and 
6  x  y 

a^  >  0  which  is  shown  in  Fig.  A- 11  requires  further  study.  The  results  of 
Sections  1  and  Z  imply  that  any  minimum  of  T  must  be  either  on  the  straight 
line 


0  <;  b  Z  <  1, 
x 

or  the  rectangle 

b  =  0, 


b 

y 


y  x 

i  -i 

x  y 


pa  I 

.  z  x 

bz  =r~=r 

X  z 


0  <  b  < 

y 


pa  I 

y  x 

i  -i 

x  y 


0  <  b 

z 


< 


I 

X 

T~ 

z 


Differentiating  the  expression  for  T  on  this  rectangle  yields  further  results 
that  all  minimum  points  on  the  rectangle  satisfy 


b 

y 


pa^ 


pa  I  b 

y  z  z 

I  +  (I  -I  )b 

y  z  y 


z 


pa  I 

„  ,  Z  X 

0  <  b  <  T - t 

z  I  —  I 


(A-  1  3) 


This  curve  is  also  shown  in  Fig.  A-ll. 


Z  Z  Z 

The  constraint  b  +  b  +  b 

x  y  z 

for  points  on  the  sphere  of  unit  radius 
In  particular,  all  minimum  points  of 


=  1  implies  that  only  values  of  T  defined 
centered  at  the  origin  are  of  interest. 

T  are  on  this  unit  sphere  so  that  for  each 
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case  of  Figs.  A-4  to  A- 11  a  minimum  of  T  can  only  occur  where  a  curve 

intersects  the  unit  sphere.  Therefore,  there  can  be  at  most  two  or  three 

minima  of  T  depending  on  the  case.  There  is  always  one  minimum  in  the 

region  b  ^  0,  b  ^  0,  b  >0  and  there  can  be  zero,  one,  or  two  minima  in 
x  y  z 

the  region  b^  <  0,  b  ^  0,  b^  >  0.  Equations  A-l  through  A- 13  may  be  used 

to  determine  the  coordinates  of  all  minima  of  T.  The  next  section  investigates 

sufficient  conditions  for  there  to  be  only  one  minimum  or  equivalently  there  to 

be  no  minima  in  the  region  b  <  0,  b  >  0,  b  >0. 

B  x  y  z 

4.  From  Eqs.  A-l  through  A-6  or  Figs.  A-l  through  A-4,  it  is  evi¬ 
dent  that  sufficient  conditions  that  no  minima  occur  in  the  region  b^  <  0, 
b  ^  0,  b  ^  0  are  as  follows. 


If  pat  =  0, 

y 

-oa  I 
x  y 


I  -I 

x  y 


>  -1. 


If  pa  =0, 
x 


pa  I 
y  x 

I  -I 

x  y 


1. 


If  pa  >  0,  pa  >0, 
x  y 


where  d 


min 


(d  .  r  >  i 

min 

is  the  minimum  value  of 


d  = 


pa  I  b 
y  x  x 

pa  I  +  (I  —  I  )b 
x  y  x  y  x 


l/Z 


the  distance  from  the  origin  to  a  point  on  the  surface  on  which  all  minima  must 

Z 

occur  in  the  region  b^  <  0,  b  ^  0,  b^  ^  0.  Differentiation  of  d  with  respect 

to  b  yields 
x  7 


(d  .  y 

min 


=  P 


,  (a  I  >2/3  ♦  (a  I  )2/3 

2  L  x  y  y  x 


(I  -I  ) 

x  y 
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In  all  three  cases,  the  sufficient  condition  for  no  minima  in  the  region  b  <0, 

x 

b  ^  0,  b  ^  0,  which  implies  that  there  is  only  one  minimum  of  T,  is 

y  z 


Z 

P  > 


(i  -i  y 

x  y 


[(a  I  )2/3  +  (a  I  )2/3l3 
L  x  y  y  x  J 


An  analogous  sufficient  condition  for  there  to  be  only  one  minimum  of  T  can 
be  derived  from  Eqs.  A-7  through  A-12.  It  is 


„2  > 


(I  -I  ) 

x  z 


fta  I  )2/3 

L  X  z 


+  (a  I  ) 

Z  X 


2/3J3 
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